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RESEARCH MEMORANDUM 

EFFECTS OF COMPRESSIBILITY ON THE FLOW PAST 
THICK AIRFOIL SECTIONS 
By Bernard N. Daley and Milton D. Humphreys 

SUMMARY 



Six, 3 -inch -chord symmetrical airfoil sections having systematic 
variations in thickness and thickness location were tested at Mach 
numbers near flight values for propeller -shank sections The tests, 
the results of which are presented in the form of schlieren photo- 
graphs of the flow past each mcdel and pressure -distribution charts 
for two of the models, were performed to illustrate the effects of 
compressibility on the flow past thick' symmetrical airfoil sections. 
Representative flow photographs indicated that at Mach numbers 
approximately 0 .05 above the critical Mach number a speed region 
was reached in which the flow oscillated rapidly and the separation 
point and the location of the shock wave were unstable. Fixing the 
transition on both surfaces of the airfoil was effective in reducing 
these rapid oscillations . The pressure distributions showed that 
the section normal -force coefficients for thick airfoils were very 
erratic at subcritical speeds; at supercritical speeds the section 
normal -force coefficients for the thick airfoils became more regular. 
Drag coefficients showed that considerable drag decreases can be 
expected by decreasing the mcdel thickness ratio . 


INTRODUCTION 


The rotational speed of the propeller has been limited hereto- 
fore by the compressibility losses at the blade tips, but as a 
consequence of the development of higher speed aircraft the thicker 
root sections having lower critical Mach numbers are now contributing 
serious adverse effects and losses more pronounced than those of 
the tip sections (reference l) . The purpose of this investigation 
was to illustrate the effects of compressibility on thick symmetrical 
airfoils suitable for use as propeller -shank sections. By utilizing 
the pressure distribution data and the schlieren photographs of the air 
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flow, force cho-racteriet-ice and information concerning the flow 
phenomena have "been detained- Test results for systematic variations 
of the maximum thickness -to-chord ratios, as well as results for 
variation of the chordwise location of maximum thickness, are included 
herein . 


SYMBOLS 


Cd section drag coefficient 

Cj section lift coefficient 

c n section normal -force coefficient 

P pressure coefficient 

P cr critical -pressure coefficient 

M Mach number 

M ch choking Mach number 

M cr critical Mach number 

q stream dynamic pressure 

a angle of attack 

* APPARATUS, MO EELS , AND TESTS 


Tunnel.- The tests were made in the Langley rectangular high- 
speed tunnel. The tunnel (fig. l) is of the closed throat, non- 
return, induction type, utilizing air at high pressure in an annular 
nozzle located downstream frem the test section to induce a flow of 
air from the atmosphere through the tunnel . The height of the test 
section is l8 inches and the width is k inches. The k-inch walls 
are flexible; adjustments of the contours of these walls produce 
variations of the longitudinal static -pressure gradient of the tunnel. 
For these tests the wall settings were such as to maintain uniform 
pressure throughout the test section without a model at any Mach 
number below 0.85. The airfoils completely spanned the k-inch 
dimension of the tunnel and were supported, in such a manner that no 
external interference to either the air flow or to the optical field 
was possible . For a more complete description of this tunnel 3 ee 
reference 2. _ _ 
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M odels . - The airfoil models investigated had 3 -inch chords 
and were as follows: 


Airfoil section 

Maximum- thi ckno s s 
location 
(percent chord) 

Leading-edge 

radius 

Reference 

MCA 0040-63 

30 

Normal 

3 . 

MCA 00 40-64 

40 

Normal 

3 

MCA 0040-65 

50 

Normal 

3 

MCA l6-040 

50 

0.444 x Normal 

4 

MCA 16-025 

50 

.444 x Normal 

4 

MCA 16-015 

50. 

.444 x Normal 

4 


The airfoil profiles are presented ir. figure 2. 

Optic al equipme nt . - The optical equipment used in those 
tests (fig. l) is similar to that described in reference 5* Changes 
in the index of refraction of the air associated with changes in 
air density are made visible with this apparatus. Compression 
shock waves and separation phenomena, as well an generalized regions 
of density change, are shown on the photographs as regions of various 
degrees of light . 

Tests ■ - Photographs of the flow around bhe alii* oil models a*t 
zero angle of attack were made by means of the schlicren method 
at all attainable supercritical speeds . Supplemented schlieren 
tests wore ms.de of the MCA 0040 -65 airfoil to determine the effects 
of fixing transition at -the quarter -chord locati on . (See figs. 3 to 10. 
High-speed motion pictures of the flow were also made by means of 
the schlieren method of air-flow photography (fig. ll) • l5ata 
obtained from these schlieren tests are presented in figure 12 . 

Pres sure -distribution data for the MCA l6-040 and the MCA 16-025 
airfoil3 were obtained at angles of attack from -2° to 12°. (See 
figs . 13 to 16 . ) Integrations of the pressure -distribution plots 
provided normal -force and pressure -drag coefficients for these 
models. The approximate Mach number range was from 0.33 to 0-73 
for the 40 -percent -thick airfoil models, from 0-33 to 0 .78 for the 
25 -percent -thick airfoil model, and from 0.33 to O.85 for the 
15 -percent -thick airfoil model. These Mach number ranges are 
comparable to those encountered in flight. The high-speed range 
was limited by the choking speed of the tunnel (reference 6) . Air- 
stream Mach number was determined by the pressure measured at a 
calibrated static -pressure orifice upstream of the test section. The 
Reynolds number range for these tests wa3 from 550 ,000 to 1,050,000. 
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PRECISION 


Constriction of the flow field of the model by the tunnel 
walls (references 6 and 7) subjects the data, to the largest source 
of error. Since constriction corrections are of questionable 
applicability at supercritical speeds, no constriction corrections 
have been applied to these data. In order to give an indication of 
the magnitude of the correction for model blockage or constriction, 
however, some corrections have been calculated by the method of 
Allen and Vincent! ( reference 7) • Although these values vary to 
some extent with model thickness ratio and Mach number, the following 
table indicates the approximate value of the correction at all sub- 
critical speeds : 


NACA 16-040 airfoil 

NACA 16-015 and 16-115 
airfoils 

Corrected c n = 0.965c n 

Corrected c = 0.975c n 

Corrected c^ - . 9650 ^ 

Corrected c^ = . 985 c ^ 

Corrected M = 1 .015M 

Corrected M - 1.010M 

Corrected q = 1.02pq 

Corrected q = 1 -015q 


Other errors due to velocity fluctuations, air-stream alinement, 
tunnel calibrations, and tunnel gradients (without model) are 
appreciably less than the constriction errors and thus for practical 
purposes are believed to be negligible. Serious humidity effects 
are not expected because all data in the present paper were obtained 
when the. atmospheric relativo humidity was below a value previously 
determined by a special investigation to be satisfactory. 

The pressure -distribution diagrams and the schlieren photographs 
are presented for Mach numbers approaching the choking Mach number 
of the turmplj the basic data for section normal -force and pressure - 
drag coefficients (figs. 17, l8, 22, and 23) are presented for Mach 
numbers up to the choking condition. The normal -force cuivefl 
(figs. 17 and l8) are dashed within 0 -04 of tho choking Mach number. 

In this range the normal -force and drag data can reasonably bo assumed 
to be affected by the choked condition. In the cross plots 
(figs. 19 to 21) data within this 0.04 range- of Mach number are 
omitted. 
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RESULTS AND DISCUSSION 


Flow Characteristics 

Schlieren photographs of the flow field of the models 
(figs. 3 to 10) were made to show the effect cn the air flow 
caused by variation of two fundamental design parameters of air- 
foils - the maximum thickness and the maximum -thickness location. 

These data indicated the probability of a rapidly oscillating flow 
in the supercritical speed range of these tests . In order to 
substantiate the existence of high-frequency oscillations occurring 
in the flow in this speed range / high -speed motion pictures of the 
flow past thick airfoils at constant Mach numbers were made by 
means of the schlieren method . A representative strip of this 
movie film (fig. ll) shows a flow of irregular nature in which the 
angle of departure of the wake, the position of the point cf sep- 
aration, and the positions of the shock waves vary at high frequency. 

From a consideration of the schlieren photographs, with the ex- 
ception of the flow about models with fixed transition (figs* 6 and 7)> 
three general flow types are encountered in the supercritical speed 
region : 

(1) A eemistable, moderately separated flew in which the 
separation points of both surfaces exhibit an ei-ratic end persistent 
dissymmetry of position over a speed range extending from some 
subcritical Mach number to a value approximately 0.05 above the 
critical Mach number M (figs. 5(a) and 8(a)) 

C A 

(2) A flow, in a higher supercritical speed range, characterized 
by relatively regular high -frequency oscillations of the rake, shock 
waves, and positions of the separation points (figs, ll end 4(b) to 4(d)) 

(3) A relatively steady separated flow originating at speeds 
greater than those at which oscillating flow ( type 2 ) is encountered 
and extending to the choking Mach number of the tunnel M^ . The 
region for type 3 flew is quite narrow for the thicker airfoil3 of 
these tests and it is, therefore, difficult to distinguish between a 
true -flow type and the influence of the choked condition of the tunnel. 
The schlieren photographs for the NACA 16-015 airfoil (figs. 10(c) to 
10(e)), however, offer good substantiation for this smooth flow at the 
higher supercritical speeds. 

Figure 12 presents for each airfoil the boundaries of the various 
types of flow encountered at supercritical speeds . The data for this 
figure were obtained from a series of individual instantaneous schlieren 
photographs taken at random Mach numbers over the speed range . The 
curves for critical Mach number and the curves bounding the oscillating- 
flow region are therefore defined only by interpolation between exact 
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teat points; however, the curves may be considered as close approxi- 
mations to the speeds at which these conditions exist. 

Variation in maximum- thickness location .- The schlieren 
photographs (figs, 3 to 5) and the graph (fig. 12) indicate that 
for airfoils having constant thickness the rearward movement of 
the location of maximum thickness from 0-3" to 0.5 -chord stations 
had the following effects: the value of the critical Mach number 

was increased from 0.51 to 0 . 58 and, consequently, the Mach number 
for shock formation was delayed by approximately the same amount; 
the Mach number for the onset of oscillating flow (type 2) was 
moved to a higher value without appreciably influencing the Mach 
number range of this type of flow (fig. 12); the violence of these 
fluctuations and the chordwise shifting of the separation point 
were reduced; the point of separation was moved rearward permitting 
a narrowing of the separated wake and the consequent indication 
of lower drag. These schlieren photographs (fig3. 3 to 5) thus 
indicate that, in addition to the increase of critical Mach number 
shown also by previous work, some slight improvement at supercritical 
speeds may be obtained by a rearward location of the point of maximum 
thickness . 

Variation of maximum thickness .- The effects of decreasing the 
maximum thickness of the 16-series sections from kO percent to 
15 percent chord were to increase the value of the critical Mach 
number from 0.57 to 0.75, to delay the occurrence of oscillating 
flow (type 2) by a similar amount, and to move the separation point 
toward the trailing edge . It is apparent also from the schlieren 
photographs in figures 8 to 10 and from the data in figure 12 that 
a reduction in the thickness of the section decreased the Mach 
number range for oscillating flow (type 2) and reduced the violence 
of the wake fluctuations . It is further evidenced by a study of the 
data for the NACA 16-015 airfoil (figs. 10(c) to 10(e) and 12) that 
the extent of the Mach number range for type 3 flow is increased with 
a decrease in maximum thickness. A decrease in the thickness ratio 
is therefore the most effective method of delaying and minimizing 
the undesirable flow characteristics at supercritical Mach numbers. 

£*lxed transition .- In an effort to determine the effect of 
fixed transition on the flow, No. 60 carborundum grains were affixed 
along the span at the quarter -chord location on the NACA 0040 -65 
airfoil. This location was chosen because photographs of the flow 
past the smooth model indicated that this location was ahead of the 
most forward location of the separation point, Schlieren photographs 
are presented for the flow with fixed transition on the upper surface 
(fig. 6) and on both surfaces (fig. 7). At the lowest speeds the 
transition stripe along the quarter-chord location seemed to produce 
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an adverse effect on the wake width and the location of the separation 
point (compare figs. 5(a), 6(a), and 7(a))* In that 3peed range 
where oscillating flew (type 2) occurred for the smooth models, the 
oscillation of the walce, the movement of the separation point, and 
the movement of the shock waves have been effectively minimized by 
fixing transition (compare figs. 5 and 7)* The flow at high Mach 
numbers appears to be improved considerably by fixing transition on 
both surfaces of the airfoil; fixing transition may produce good 
results in cases in which vibrations are a serious problem* At high 
Reynolds numbers or with airfoils having the transition point already 
far forward, the oscillations of the flow may be less severe than 
found in these experiments . 


Airfoil Force Characteristics 

Airfoil normal forces .- Figures 13 to l6 show variations of 
representative pressure distributions with Mach number for the 
NACA 16-025 and the NACA 16-040 airfoils. The data of figures l4(a) , 

14(b) and 16(a) to 16(d) indicate a pronounced variation with Mach 
number of the surface loading near the trailing edge of the model. 

This variation is caused at subcritical speeds by the semistable 
unpredictable nature cf the separation point characteristic of 
type 1 flow. (The critical -pressure coefficient P cr is represented 

by the horizontal broken lines on the pressure -distribution charts.) 

These data indicate correspondingly large and unpredictable changes 
in section normal -force coefficients with Mach number. Such indications 
of changes in normal -force characteristics are confirmed by corres- 
ponding integrated section normal -force data of figures 17 and l8. 

The normal -force coefficients at subcritical and low supercritical 
speeds where erratic separated flow (type l) predominates are generally 
irregular and these data, probably cannot be reproduced accurately . 

Also as a result of assymetrical separation phenomena, the angle of 
attack for zero lift is considerably displaced . A similarly irregular 
zero -lift angle has been observed in data for the NACA J.6-539 airfoil 
as tested in the Langley 24 -inch high-speed tunnel (reference 4). 

Those unpredictable semistable normal -force characteristics (type 1 flow) 
exhibited by the thicker sections make the application of these air- 
foils to propeller -shank sections extremely undesirable and uncertain 
in this speed range , 

As mentioned in the section entitled ’’Flow Characteristics," at 
Mach numbers approximately 0.05 above the critical, the flow becomes 
oscillating (type 2) and the separation point loses the former- 
erratic and semistable nature of type 1 flow. This speed range 
corresponds closely to the speed at which the normal -force curves (figs. 17 
and l8) first assume a regular variation with Mach number or angle of 
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attack. In this Mach number range the angle of attack for zero 
lift approaches zero. All pressures that are recorded when the 
flow is oscillating rapidly arc average pressures, and therefore 
the large flow variations that appear in the schlieren photographs 
in this speed range do not have a corresponding effect on the 
pressure data. The aforementioned phenomena are also clearly 
illustrated in figures 19 and 20, which present the variation of 
section normal -force coefficient with angle of attack for several 
Mach numbers . Lift data for the NACA 16-115 airfoil from the 
Langley 24 -inch high-speed tunnel are presented in figure 21 for 
the purpose of illustrating the improvements that can be obtained 
by using thinner airfoil sections. A general improvement can be 
observed in the regularity of the lift -curve slope The angle of 
attack for zero lift also does not vary in the Mach number range 
presented. 

Airfoil drag characteristics . - The lack of pressure recovery 
at the higher speeds shown in the pressure distributions is an 
indication of separation of the flow and, consequently, large 
pressure drags are to be expected. Pressure -drag coefficients for the 
NACA 16-025 and the NACA 16-040 airfoils were obtained by combining 
the drag components of the normal -force coefficient with those of 
the chord -force coefficient. The results are presented in figure 22 
together with total -drag data on a 5 - inch-chord NACA J.6-115 airfoil 
from tests in the Langley 24 -inch high-speed tunnel. The drag 
coefficient for the thin model is small and constant throughout most 
of the Mach number range . The pressure -drag coefficients for the 
thicker models are considerably higher in the low Mach number range 
because of severe separation effects ; these coefficients increase 
rapidly at the higher Mach numbers because of the compression shock 
losses and increased separation losses which occur as a result of 
the lower critical Mach nurnoer of the thicker sections. (See figs. 3 
to 10 and 22.) . A comparison of the data for the ceveral angles of 
attack of the NACA 16-025 airfoil shows that the drag coefficient 
for this model is comparatively insensitive to angle of attack in 
the lower Mach number range; at higher speeds the drag is markedly 
affected by angle -of -attack variations. In the Mach number range 
where a large force break i3 observed in the section normal -force 
curve (fig. 18) for the NACA l6-040 airfoil at -2°, the drag curve 
(fig 22) has a similar jog. This jog can be attributed to the 
peculiar separation phenomena which influenced the normal -force 
characteristics . 

Drag results were not obtained to show the effect of the variation 
of the location of the maximum thickness, because these effects are 
known to be of less importance than the effects of the large thickness 
changes studied . Drag data are presented primarily to show the effect 
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of variations in the maximum thickness of the airfoil. In order to 
illustrate these drag variations better, the data have been recomputed 
on the basis of frontal area. These data (fig. 23) illustrate that 
the drag of the NACA l6-040 airfoil can be decreased more than 
90 percent at a Mach number of 0 .69 by using an NACA 16-015 airfoil 
of the same maximum thickness . It thus appears that the best method 
of reducing the propeller -shank drag is to reduce the thickness-to- 
chord ratio of the propeller -shank section to a value insuring favorable 
flow characteristics . 


CONCLUSIONS 


Tests to determine the effects of compressibility on the air 
flow and force characteristics of thick airfoil sections in a Mach 
number range comparable to a range of flight values for propeller- 
shank sections indicated that : 

1. At Mach numbers approximately 0.05 above the critical value 
for thick airfoils, a speed region was reached in which the flow 
oscillated rapidly. In cases in which serious flow oscillations 
occurred, the flow was improved considerably by fixing transition 
on both surfaces of the airfoil. 

2 . Improvement in section characteristics may be expected from 
moving the maximum-thickness location from 0.3 to 0 .5 chord due to 
drag decreases caused by the increase in critical Mach number and 

by the narrowing of the separated wake, particularly at supercritical 
speeds . 

3. Relatively large improvement in section characteristics 
resulted from decreasing the maximum thickness of the 16 -series 
sections from 40 percent chord to 15 percent chord. The pressure - 
drag coefficient based on model thickness was reduced by as much as 
90 percent at a Mach number of 0 , 69 • 


Langley Memorial Aeronautical Laboratory 
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Langley Field, Va. 
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FIGURE I.- SCHEMATIC DIAGRAM OF LANGLEY RECTANGULAR 
HIGH-SPEED TUNNEL AND SCHLIEREN APPARATUS. 
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Figure Z .— Pro fit es of oirfoits Jested in the Long/eg 

rectangular h igh- speed funnel- 
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Fig. 3 



(a) M = 0.564. 


(b) M = 0.585. 



(c) M = 0.609. (d) M = 0.706. 


Figure 3.- Schlieren photographs of flow over NACA 0040-63 
airfoil. M cr = 0.510; M c ^ = 0.714; a = 0°. 
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Fig. 4 


(b) M = 0.626. 




Figure 4.- Schlieren photographs of flow over NACA 0040-64 airfoil. 

M cr = 0.550; M ch = 0.721; a =0°. 


N AT 10 N AL ADVISORY COMMITTER FOR AERONAUTICS 
LANGLBT MEMORIAL AERONAUTICAL LABORATORY - LANGLEY PI ELD, VA. 


NACA RM No. L6J17a 


Fig. 5 



(a) M = 0.602. 


(b) M = 0.678. 




(c) M = 0.707. 


(d) M = 0.731. 


Figure 5.- Schlieren photographs of flow over NACA 0040-65 airfoil. 

M cr = 0.580; Mc h = 0.731; a =0°. 
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Fig. 6 




(a) M = 0.600. 


(b) M = 0.677. 



(c) M = 0.707. (d) M = 0.729. 

Figure 6.- Schlieren photographs of flow over NACA 0040-65 airfoil 
with transition fixed at 0.25 chord on upper surface. M cr z 0.580; 

M ch = °* 729 ’ a = °°* 
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Fig. 7 




(a) M = 0.600. 


(b) M = 0.676. 



Figure 7.- Schlieren photographs of flow over NACA 0040-65 airfoil 
with transition fixed at 0.25 chord on the upper and lower surfaces 
M cr = 0.580; M ch = 0.727; a = 0°. 
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Fig. 8 




(c) M = 0.711. 


(d) M = 0.723. 


Figure 8.- Schlieren photographs of flow over NACA 16-040 airfoil. 

M cr = 0.575; M ch = 0.729; a = 0°. 
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Fig. 9 



(a) M = 0.674 . 


(b) M = 0.706 . 




(c) M = 0.758. 


(d) M = 0.776. 


Figure 9.- Schlieren photographs of flow over NACA 16-025 airfoil. 

M cr = 0.657; M ch = 0.787; a = 0°. 
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Fig. 10 




(b) M = 0.788. 


(a) M = 0.756. 




(c) 



(d) M = 0.826. 


(e) M = 0.838 . 


Figure 10.- Schlieren photographs of flow over NACA 16-015 airfoil. 

M cr = 0.730; M ch = 0.847; a = 0°. 
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Fig. 11 



Figure 11.- Continuous motion -picture sequence illustrating 
oscillating flow. Camera speed, 2060 frames per second; 
NACA 0040-63 airfoil; M = 0.63; M cr = 0.510; M ch = 0.714; 
a = 0°. 
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Fig. 12 



Figure I2r Supercr i t'ical flow bound a ri es 
for NHCH airfoils tested. oc=0°. 
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P/gure iJ. — Pressure distri bulion for an NACA /6'025 
Q/rfo/Z. /% r =o. 6 S 7 ; o. 787 ; a=o°. 
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Fig. 14 
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P/gure /4 . — Pressure distribution for an NACA 16-025 
Q/rfo/Z. P7 cr = O. L4L ; OJ&5 ; CC s 4° . 
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Fig ure 15. — Pressure distribution for an NACA \6~0A0 
a/rfo/Z. AZ c/ u O.S75 ; =0.7Z? ; OC= 0\ 
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Fig. 16 
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Pi pur e id . — Pre ssure distribution for an NACA 16 040 
a/rfo/Z. /d^ r =O.S44 ; /Z^p^o.787 4 oc- -a e . 


Fig. 17 
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Figure 


17.- Variation of section normal-force coefficient 
with Mach number. NACA /6-OZ5 airfoil. 
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/Soure Hr Variation of section normal-force 
coefficient with anqie of attacK. 
NACA / 6-025 airfoil. 
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Fig. 20 
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Figure 20 — Variation of section normal- force 
coefficient with anqie of atfacr. 
A/ACA /6-040 airfoil . 
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